On the Limb Darkening of Planetary 
Atmospheres in the Thermal Infrared 


Thermal maps of Venus (Murray, et al . , 1963; Westphal , et al., 1965) 
and Jupiter (Murray, et al., 1964) have demonstrated that limb darkening 
in the 8-14 micron wavelength interval is pronounced for these two 
planets. Since both these planets are completely cloud-covered, it 
is possible that both anisotropic multiple scattering (Goody, 1964) and 
a temperature gradient (Pollack and Sagan, 1965) are responsible. It 
will be the purpose of this note to demonstrate how exact solutions to 
the standard problem- -that of describing the angular distribution of 
outgoing radiation from a slab of material bounded on both sides by a 
vacuum (Chandrasekhar, 1950) --may be expressed in terms of Chandrasekhar's 
•if- and 0- functions of zero order when both scattering and temperature 
gradient effects are present. The additional problem of imposing the 
correct lower boundary conditions to account for surface effects in 
order to specify the angular distribution of outgoing thermal radiation 
(limb darkening) from the atmosphere is quite straight-forward (cf. 
Chandrasekhar, 1950, pp. 269-274) and will not be presented here. 

Let Ip ( t , p.) be the monochromatic specific intensity of thermal 
radiation in a doubly infinite plane-parallel atmosphere. Consider 
an atmosphere of finite thickness t x = t 3 - t 3 having scattering and 
thermal emission properties identical with those th - the doubly infinite 
atmosphere has in the range (t 2 £ t s t 3 ). Let l(0,+n) and X ( t-j * — jju J 
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designate the outgoing radiation fields in the directions +p, and -p, 

(0 < p, ^ 1) at the levels x = 0 and x = t, respectively. 

By developing an invariance principle known to Chandrasekhar (1950), 
Mullikin (1962) and many others, it is found that the solutions for 
the standard problem are given by 
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Hx, r/ u) = Z f <VW -±f S(r l}ri y)I p (r, iy ,')^ 
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where 

1(0, -p,) = I( Ti ,+|J.) = 0 (0 < n s 1), 

and where S and T are the diffuse scattering and transmission functions 
in Chandrasekhar's (1950) notation for atmospheres of normal optical 
thickness x x • Physically (1) and (2) arise from the fact that any 
mass element at a level x(0 s x ^ h) is unable to distinguish between 
being imbedded in a doubly infinite atmosphere, and being imbedded in a 
finite atmosphere that is irradiated on the top and bottom by the radia- 
tion fields l p (0,-p,) and I p (x 1 , + u) respectively. 
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Consider the equation of transfer for a partially thermally 
emitting, partically anisotropically scattering plane-parallel atmos- 
phere: 

+/ 

~ i f P I ( r j/) ft ' - (l-tto)Bfr-) . (4) 
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Let pOiSli.p- 1 ) = p(g,,p.') be independent of j* Further let p ( jj, , ^ 1 ) 
be the azimuth-independent part of the expansion in spherical harmonics 
of a finite Legendre polynomial series approximation of the phase func- 
tion for single scattering, and let B( T ) be approximated by a finite 
power series in x; i.e., 
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It can be shown (cf. Samuelson, 1965) that a particular integral 
satisfying equation (4) under conditions (5) - ( 6 ) is given by 
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where the constant coefficients c are uniquely generated by the 
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That equation (7) is also the complete solution for a doubly 
infinite atmosphere is rendered plausible on physical grounds by noting 
that [cf . eq. (6) ] 
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i.e. the limiting values of the average specific intensity and the net 
flux are respectively equal to and proportional to the limiting values 
of the Planck function and its divergence. 


Now Chandrasekhar's i|r- and 0- functions of order zero and degree 

are defined by (Chandrasekhar, 1950) 

x J J . . 
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Upon multiplying through both numerator and denominator in all the 
integrands of equations (1) - (2) by p, 1 , it is readily demonstrated, 
with the aid of equations (7) and (12) - (13) and the relation 

that (1) and (2) respectively reducfe to 
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where all i|r and 0 of negative degree are defined to be bounded in the 
interval (0 £ p, £ 1). 

Equations (15) - (16) are notable in two respects. First, the 
only non-elementary functions involved are the \|r- and 0- functions of 
zero order. For a given and set of = 0,..,N[cf. also eq. ( 9) ] J" 

these functions may be calculated once and for all with the aid of existing 
computer programs (cf. Churchill, et al., 1961). And secondly, the 
complete description of the thermal properties of the atmosphere are 
contained entirely in the constant coefficients c„ (r = 0,...M; s =' r,...M). 
In view of the foregoing remarks and the simplicity of relations (8), 
the practical problem of solving (15) - (16) numerically should not be 
inordinately difficult. The additional problem of calculating the planetary 
limb darkening under restrictions (5) - (6) in accordance with Chandrasekhar's 
discussion of the "planetary problem" (loc. cit.) is quite straight- 
forward. 


R. E. Samuel son 

Goddard Space Flight Center 

Greenbelt, Maryland 
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